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ABSTRACT 

We apply a simple statistical method (Derenzo & Hildebrand 1969) to 
estimating the completeness of quasar surveys. It requires that an area has been 
covered by two or more, preferably different, selection techniques. We use three 
suitable data sets with separate selections from: variability and UV-excess (170 
quasars); objective prism and UV-excess (141 quasars); multicolour and X-ray 
(ROSAT, 19 quasars). We find that, for selection by UV-excess, the common 
limit of U — B < —0.35 ± —0.05 leads to losses of ~ 35%, typically missing 
low-luminosity {Mb ^ —24.5) quasars, independently of redshift. Systematic 
incompleteness will therefore affect the new generation of large quasar surveys 
that select by U — B < —0.35. By correcting for this incompleteness, we find, 
from the first data set (B < 21.0 and z < 2.2), that the evolution of the quasar 
luminosity function (LF) is best described by joint luminosity and density 
evolution. When extrapolated to z — 0, the LF matches that of local Seyfert 
galaxies better than any previous determination. The LF shows an increase in 
the number of low-luminosity quasars at low redshifts and of brighter quasars 
at intermediate redshifts, relative to the LF of Boyle et al. (1990). This result 
is consistent with models in which quasars fade from an initial bright phase. 

Subject headings: surveys — quasars: luminosity function — galaxies: Seyfert 
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1. INTRODUCTION 

The completeness of a quasar survey should properly be specified by the survey 
selection function, P(M,z,SED), which gives the probability of detecting 
function of absolute magnitude, redshift and spectral energy distribution (SED). The 
number of detected quasars in each cell of (M, z, SED) space is divided by P(M, z, SED) 
to give the number of quasars that could be found. Summation over all cells gives the total 
population of quasars within the survey limits. In practice, P(M, z, SED) can not usually 
be determined reliably, and many cells may be empty or too poorly occupied for reliable 
estimates of the true number. 

In this paper, we apply a statistical method to estimating the completeness of quasar 
surveys. It requires that an area of sky has been covered by two or more, preferably 
different, selection techniques. It does not require P(M, z, SED). 

We adopt H = 50 km s" 1 Mpc -1 , q = 0.5 and a = 0.5 (f u oc v~ a ). 



2. ESTIMATING COMPLETENESS: THE METHOD 

The method for estimating the completeness of a quasar survey involves the comparison 
of the number of detected quasars and an estimate of the total population. The total 
population is estimated using the method devised by Perenzo fc Hildcbrand (1969) for 



estimating scanning efficiencies in the detection of events in particle physics. The use of 
the Derenzo & Hildebrand (subsequently DH) method here is appropriate because in both 
applications the size of a total set of events is estimated from repeated scans of the set, with 
each scan yielding only a subset. In DH, the total number of particle events is estimated 
from repeated scans of film records; in this paper, the total number of quasar events is 
estimated from repeated scans (surveys) of an area of sky. 
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In a similar application to ours, Harwit & Hildebrand (1986) used the DH method 
to estimate the remaining number of observational phenomena still to be discovered 
in astronomy. Four kinds of scans were defined according to their usage of (i) optical 
continuum methods; (ii) radio techniques; (iii) optical spectroscopy; and (iv) all other 
techniques to discover astronomical phenomena. 

As a demonstration, we have applied the method to two samples of objects derived 
from a known parent population. We show that the total number of objects in the parent 
population is satisfactorily recovered. 

The net visibility, v, of an object is defined to be the probability that it will be found by 
an average detection technique. It takes values in the range (invisible) to 1 (unmissable). 
Note that the net visibility depends on both the intrinsic characteristics of the object and 
the percept iveness of the detection techniques. 

The net visibility function, F(v), is defined such that F(v)dv is the number of objects 
with net visibilities in the range v to v + dv. The total number of objects is then 



Note that F(v) is an average of visibility functions for the individual techniques 
(Derenzo & Hildebrand 1969). The existence of visibilities and visibility functions is an 
assumption of the method. See the Appendix for a formal derivation of F(v) and the above 
result for N T . 

If all objects had equal prior probability of being detected then F{v) would be a 
5-function at the average detection probability. In general, however, F(v) is an unknown 
function that must be estimated. 

Consider a population that is surveyed using several techniques, preferably different 
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but not necessarily. Suppose, for example, that there are three techniques. A first technique 
(A) detects A objects, another technique (B) detects B objects, and a third technique (C) 
detects C objects. A particular object may be detected by more than one technique. 

The average number of objects detected by one application of the average technique, 
M[l], is the average number detected by the three techniques 



To determine the average number of new objects detected by the second application of the 
average technique, we take the techniques in pairs. A new object is detected if it is seen by 
one technique but not by the other. Thus, 



where Ab denotes the number of objects detected by technique A but missed by technique 
B, and so forth. Similarly, the average number of new objects found by the third application 
of the average technique is 



In general, the distributions of undetected objects after one, two and % — 1 surveys 
are respectively F(v )(1 — v), F(v)(l — v) 2 and F(v)(l — v) 1-1 . Hence the number of new 
objects found in the i th survey using the average technique will be 



The collection of values M[l], M[2], . . . , M[n] can then be used to evaluate F(v ), if we 



M[l] = [A + B + C]/3 



M [2] = [Ab + Ac + Ba + Bc + Ca + Cb]/6 



M[3] = [Abc + Bac+Cab]/3 




(1) 
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assume for F(v) an expression with not more than n parameters, where n is the number of 
surveys. 



Harwit & Hilde brand (1986) argue that, in general, the total number of objects, Nt, 
derived from F(v) will be approximately the same for any smooth trial function, F(v), 
whose m parameters (m < n) are adjusted to give satisfactory fits to the values of M[i]. 
Derenzo fc Hildebrand (1969)| show that for any function F(v) that is integrable in the 
region < v < 1, and for all values of M[i], i = 1, . . . , oo, which satisfy equation [[], a lower 
limit for Nt is given by 



N T > 



Em 



+ 



M 2 [n] 



M[n] 



i -i 



M[n - 1]_ 



subject only to statistical fluctuations in the measured quantities M[l], . . . , M[n]. 
In the case where n = 2, this expression can be rewritten as 



N T > 



;aub + ahb) 2 

4(A fl B) 



This shows that, when there are only two surveying techniques, the estimated lower limit 
for Nt depends on the number of objects that both techniques have in common (A fl B). 
When the total number of objects detected by both techniques is very large but the number 
of objects in common is relatively small, this dependency will lead to a very large estimate 
of Nt- From an analysis of the behaviour of the above expression in these circumstances, 
we have found that unreliable estimates of Nt tend to be given when (A U B)/(A D B) ^ 5. 
We therefore recommend that the method is applied only to data sets in which this ratio is 
<^ 5. There are no such discontinuities for higher values of n. 

As a demonstration of the method, we consider the optical survey by Tritton & Morton 
(1984). They identified all objects in a region of 0.31 deg 2 , finding 747 ordinary stars, 3 
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white dwarfs, 4 quasars, and 143 galaxies. We take the 601 stars found with B < 20.0 
as a complete sample and then isolate two subsets by different selection criteria: (i) 379 
stars with B — V < 1.0; (ii) 312 stars with an objective-prism classification of spectral type 
earlier than K. Subset 1 has 111 stars not found in subset 2; subset 2 has 44 stars not 
found in subset 1. This gives M[l] = 345.5 and M[2] = 77.5 and hence > 445. The 
trial functions used in the following sections give estimates of the total number in the range 
445-667, which embraces the actual total of 601 stars. The median estimate is 546. 



3. ESTIMATING COMPLETENESS: THE SURVEYS 

The method estimates the total number of quasars in the parameter space defined 
by the boundaries of a data set — that is, the total number of quasars lying within the 
magnitude, redshift and SED ranges of the data set. We have applied the method to 
three suitable data sets, each based on two different survey techniques but satisfying our 
condition that (A U B)/(A n B) 5. 



Data set 1 is based on a survey of quasars in ESO/SERC field 287 by [Hawkins fc Veron 



(1995)| (see also: Hawkins fc Veron 1993| ; Veron fc Hawkins 1995| ). Quasar candidates were 



selected primarily by variability (s > 3.5, where s is a variability parameter typically ten 
times the amplitude of the variability), although a UV-excess (UVX) criterion (U — B < 0) 
was also used. The survey actually consists of three samples: (i) a faint sample of 111 
quasars covering 4.6 deg 2 with 19.5 < B < 21.0 and z < 2.2; (ii) a bright sample of 117 
quasars covering 18.8 deg 2 with B < 19.5 and z < 2.2; (iii) and a high-redshift sample of 86 
quasars covering 18.8 deg 2 with B < 21.0 and z > 2.2. We construct data set 1 from the 
faint sample and those quasars from the bright sample that lie within the boundaries of the 
faint sample. It comprises 170 quasars covering 4.6 deg 2 with B < 21.0 and z < 2.2. We 
have not considered quasars with z > 2.2 as they were not sampled by two independent 
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methods but were selected jointly on their colour and variability. The intrinsic luminosity 
limit for data set 1 is Mb = —19.83. 145 quasars would be detected by their variability 
(s > 3.5) and 163 by their colour (U — B < 0), giving M[l] = 154. Only 7 quasars are found 
by variability that are not found by UVX whereas UVX finds 25 quasars not detected by 
their variability, which gives M[2] = 16. 

Data set 2 is based on two quasar surveys of ESO/SERC field 927: (i) a sparse-sampling 
survey of 118 quasars, selected by the Automated Quasar Detection software (AQD; |Clowe; 



19(8|; |Clowes, Cooke fc BeardT984l ), covering 25deg 2 with B < 20.8 and z < 2.8 (|Clowes| 



fe |Campusano 1994] ; praham 1997| ; Clowes, Campusano & Graham, in preparation); (ii) 
a survey of 151 quasars, selected by UV-excess U — B < —0.3, covering 14.4 deg 2 with 
B < 20.0 and z < 2.8 flGraham 1997| ). Data set 2 is constructed from all the UVX-selected 
quasars with z < 2.2 and all the AQD-selected quasars with B < 20, z < 2.2 that lie 
within the boundaries of the UVX survey. It comprises 141 quasars in 14.4 deg 2 and has an 
intrinsic luminosity limit of Mb = —21.04. 135 quasars were found in the UVX survey and 
49 quasars in the AQD survey satisfying the relevant criteria. Only 6 quasars were found 
by AQD that were not detected by UVX whereas UVX found 92 quasars not detected by 
AQD. Thus, M[l] = 92 and M [2] = 49. 

Data set 3 is based on two surveys of Selected Area 57 (SA57): (i) a survey of 30 
quasars, selected by multicolour techniques, covering 0.29 deg 2 with B < 22.6 and z < 3.12 
( 1-K-oo fc Kron 1958] ); (ii) a survey of 19 quasars, selected from ROSAT HRI observations 



of SA57, covering 0.13 deg 2 with a limiting flux of 0.2 x 10 14 erg s 1 cm 2 and z < 3.02 
( |Miyaji et al. 1997| ). Data set 3 is constructed from all .ROS^T-selected quasars and all 



quasars from Koo & Kron that lie within 12' of the ROSAT pointing centre. It comprises 
19 quasars covering 0.13 deg 2 . 17 quasars were identified by Koo & Kron that meet the 
relevant criteria and 19 by ROSAT. No quasars were found by Koo & Kron that were 
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not found by RO SAT and 2 were found by RO '5 AT 'that were not found by Koo & Kron. 
Thus, M[l] = 18 and M[2] = 1. Note, however, that of the 33 sources detected by ROSAT 
in SA57, 7 are classed as faint stellar objects. These could be quasars but they would be 
fainter than the magnitude limit of Bj = 22.5 of Koo & Kron flMiyaji 19971 ). 



We have used the values of M[l] and M [2] for the three data sets to fit different trial 
functions, F(v). The F(v) and the corresponding estimates of the total number of quasars 
are given in Table [l]. The figure in parentheses after each estimate is the completeness of 
the data set based on the estimated total number for the data set. The lower limit to the 
total number of quasars is also given. With two techniques contributing to a data set the 
trial functions have only two free parameters. However, the trial functions give a good 
coverage of the possible functional forms of the true visibility function. 

The estimates show that: data set 1 is 63-99% complete; data set 2 is 47-72% complete; 
and data set 3 is 61-100% complete. If we use the median estimates then: data set 1 is 89% 
complete; data set 2 is 54% complete; and data set 3 is 90% complete. If we assume that 
there were errors of \/N in the numbers of common (A, B) and distinct (Ab, Ba) objects 
comprising each survey then the typical error in the estimate of the total number of objects 
is ~5% in data set 1, ~5% in data set 2, and ~15% in data set 3. 

The success rates, for median estimates, of the individual selection techniques are: 
variability (with s > 3.5) selects 76% and UVX (with U - B < 0) 86% of the total 
population of quasars for data set 1; UVX (with U — B < —0.3) selects 52% and AQD 19% 
of the total for data set 2; (the AQD survey is by sparse sampling); colour selection finds 
81% and X-ray 90% of the total for data set 3. 
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4. COMMENTS ABOUT SURVEYS 

Random incompleteness — effectively sparse sampling — does not affect the 
expectation values of any statistics if the sampling rate is known ( |Kaiser 1986Q . However, 



systematic incompleteness (e.g. bias against the detection of low-redshift or low-luminosity 
quasars) means that a sample is non-random and incorrect conclusions may be drawn. 

Samples of quasars selected by UV-excess (U — B < —0.35 ± 0.05) are normally 
assumed, following [Veron (1983)| , to contain at least 90% of the quasars with z < 2.2 



within a given survey. However our results suggest that a large fraction of quasars with 
z < 2.2, possibly as high as 50% (from data set 2), do not have a large UV-excess and must 
have been missed by previous UVX surveys, e.g. [Boyle et al. (199U}| . Such losses could 



have serious consequences if the incompleteness is systematic. We can test for systematic 
incompleteness by examining the effectiveness of selection by UVX for different U — B 
limits. To do this, we take data set 1 and consider UVX limits of U — B < —0.35 and 
U-B < 0.0. 

One possibility is that a limit of U — B < —0.35 is biased against low redshift quasars 
{z < 0.7). A plot of U — B against redshift for all quasars in data set 1 suggests that this is 
likely (see Hawkins & Veron 1995, Fig. 1). For data set 1, Table p] shows the effectiveness 
of selection by UVX for z > 0.7. The figure in parentheses after each estimate is the 
effectiveness of selection by UVX with the relevant UVX limit. The lower limit to the total 
number of quasars present is also given. If only low redshift quasars were being missed 
for U — B < —0.35 then for z > 0.7 the effectiveness for the two UVX limits should be 
approximately the same. This is clearly not the case. 

A second possibility is that a limit of U — B < —0.35 is biased against low 
luminosity quasars. Several authors have considered such a bias as likely but estimate the 
incompleteness at ~ 10% (e.g. Boyle et al. 1990). Note that a survey that arbitrarily defines 
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a quasar to have Mb < —23 (Schmidt & Green 1983) would also be biased against low 
luminosity quasars. For data set 1, Table [I] similarly shows the effectiveness of selection by 
UVX for Mb < —24.5. If only low luminosity quasars were being missed for U — B < —0.35 
then for Mb < —24.5 the effectiveness for the two UVX limits should be approximately the 
same. This is indeed the case. 

In summary, a quasar sample selected by U — B < —0.35 is systematically incomplete, 
being biased in particular against low luminosity quasars. Part of the incompleteness is 
probably caused by detectable radiation from the host galaxy. This will redden the quasar 
and it will cause losses in surveys that require candidates to be classified as star-like. 



From Parkes flat-spectrum quasars, [Francis et al. (1997)| have suggested the existence 



of three categories: (i) conventional blue quasars (40%); (ii) red quasars, with reddening 
from synchrotron radiation (30%); (iii) very red quasars, with reddening from dust (30%). 
These last two categories of red quasars, being faint in B, could be the type of objects that 
comprise the low luminosity population suggested by our analysis. 

The incompleteness revealed by this work, has implications for some of the new 
generation of large quasar surveys. For illustration, consider the AAT 2dF quasar survey 
flSmith et al. 1996|) . This currently specifies U -B < -0.36 and 18.25 < B < 21.0, and has 
a projected final completeness of ~ 90%. Adapting data set 1 to have the same selection 
criteria gives M[l] = 124 and M[2] = 29.5. The corresponding lower limit to the total 
number of quasars is > 163 and the median value of the five estimators is 217. 104 
quasars are detected with U — B < —0.36. Thus, our results imply that at best the AAT 
2dF survey will select only 64% of all quasars with 18.25 < B < 21.0 and z < 2.2. 
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5. 



THE LUMINOSITY FUNCTION REVISITED 



Previous determinations of the quasar luminosity function (LF) could have significantly 
underestimated the density of low-luminosity quasars. We can now estimate the 
incompleteness, correct for it, and so determine the 'true' LF. 

For data set 1, we have divided the ranges of apparent magnitude and redshift 
into bins: 16 < B < 18, 18 < B < 19, 19 < B < 20, 20 < B < 21; < z < 0.7, 
0.7 < z < 1.5, 1.5 < z < 2.2. Within each joint magnitude and redshift bin, we have 
estimated the 'true' surface density of quasars using both the maximum completeness 
estimator, F(v ) = K5(v — a), and the median completeness estimator, F(v) = K(l — v) a . 
We have generated 'best-fit' models to these data using maximum likelihood techniques for 
various parametric forms for the LF. We consider two models involving pure luminosity 
evolution with parametric forms given by Poyle, Shanks fc Peterson (1988)] : a smoothed 
double power-law function 



where, in both cases, the redshift dependence, M(z), can be expressed either as a (1 + z) 
power-law evolution, M(z) = M* — 2.5fcilog(l + z), or as an exponential evolution with 
look-back time, M(z) = M* — 1.08/clt. The look-back time expressed as a fraction of the 
age of the universe is r(z) = 1 — (1 + z)~ 3 / 2 . $* is a normalisation factor. 

We consider also a model with density and luminosity evolution parameterized as 



$(M,z) 



IQ0A(M-M(z))(a+1) _|_ IQ0A(M-M(z))(f3+1) 



and a single power-law function 



$(M,*) 



]_Q0.4(M-M(z))(q+1) 



$(M,z) 



IQ0A(M-M(z)){a+1) + IQQA{M-M(z))(0+1) 
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where the density evolution function is expressed as an exponential with look-back time 

p D (z) = exp(k D r) 

a, (3, M*, ki and ko are the free parameters that determine the fit. 

The maximum-likelihood solution is found in each case by minimizing the function: 

S = -2 x ij lo S e + 2 Ikj 

where is the observed surface density in magnitude bin i and redshift bin j, is the 
expected surface density in bin i and bin j, and the summation is over all bins. 

Each 'best-fit' model was tested for goodness-of-fit using the Pearson x 2 statistic. 
Table [2] gives the best-fit parameters for the models and their x 2 probabilities, p(> x 2 )- For 
comparison, we have calculated how well some of the best-fit models of [Boyle et al. (199U)| . 
normalised appropriately, fit the data; these are also included in Table |2|. 

We can reject a single power-law LF and all of the best-fit models of Boyle et al. Of 
the remaining models, all of which incorporate double power-law LFs, the data are well 
described by those with density evolution. The best fits are obtained when all evolution is 
expressed by exponentials with look-back time. Fig. 1 shows our best-fit models applied 
to the data; for comparison, the best-fit model of Boyle et al. (model B), normalised 
appropriately, is also shown. At low redshifts the Boyle et al. model underestimates the 
number of faint quasars and at intermediate redshifts it underestimates the number of 
brighter quasars. Fig. 2 shows the LF corresponding to our best-fit model extrapolated to 
z = 0. For comparison the figure also shows: the LF derived from local Seyfert galaxies of 
types 1 and 1.5 (Cheng et al. 1986); the LF for soft X-ray selected AGN ( |Franceschini et aT. 
19(93]), extrapolated to z = 0; and Boyle et al. model B, extrapolated to z = 0. Our model 
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provides the best-fit to the local Seyferts. The curve for X-ray selected AGNs (including 
all types) lies a factor of ~ 2 above our LF, which suggests that other types of AGN (e.g. 
starburst galaxies) are not related to quasars. 

Compared with the most successful model of Boyle et al. (model B), our LF shows an 
increase in the number of low luminosity quasars at low redshifts and of brighter quasars at 
intermediate redshifts. It also matches well the LF of local Seyfert galaxies. Such behaviour 
is consistent with quasar models in which evolution is caused by the progressive exhaustion 
of the fuel supply to the central black hole. Quasars fade from an initial bright phase until 
a low, quasi-steady rate of energy production is reached. This rate declines only slowly over 
a long timescale. Quasars in this final phase of evolution are equivalent to local Seyfert 
galaxies. 

The referee provided helpful comments. MJG acknowledges the support and hospitality 
of Departimento d'Astronomia i Astrofisica at the Universitat de Valencia. We thank Dr T. 
Miyaji for providing us with data before their publication. LEC was partially supported by 
FONDECYT grant number 1970735. 

APPENDIX 

In this appendix, we show that the total number of objects can be obtained from an 
average visibility function. 

Consider a data set that is scanned N times. Each scan has its own visibility function, 
fiiv). Let Xi denote the number of objects found by scan i and Xi be the number of objects 
not found by scan i then, for example, the expression XiXjXk represents the number of 
objects found by scan i but not by scan j or scan k. 
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Consider the average number of objects found in one scan, 

1 N 

i 

In terms of the visibility function for a scanner, 

/ vfi(v)dv 
Jo 



X, 



and so 

L N r l 



1 r 1 N 

If we define the average visibility function by 



1 N 



N 

then 



M[l] = C vF(v)di 
Jo 



Now consider the average number of new objects found in the second scan, 



i 



v ' 1,3(1^3) 



l x j 



For a single scanner (fixed i), 

1 N f 1 
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and so 

l N r 1 



M P] = ^E/ q vf t (v)(l-v)di 



1 r 1 N 



NJo 



[ v^Mv)(l-v)dv 



= f 1 vF(v)(l-v)di 
Jo 



In general, the average number of new objects found in the i th scan will be 
M[i] = 

T,^j,...,k(i^...^k) x i x j ■ ■ ■ X k 



N(N-l)...(N-i+l) 



For a single scanner (fixed i), 



1 



(N — 1) . . . (N — i + 1) , „. , XlXj 



j,...,k(j^-¥k^i) 

° l vf % {v){l-vY~ l dv 







and so 



M[i] = C vF(v)(l -vY^dv 
Jo 

Now consider the integral, 

r i 1 r 1 N 

I Hv)dv = £/,<„)* 

1 N /"I 



Thus, fitting an average visibility function to the data allows the total number of objects to 
be estimated. 
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Fig. 1. — This shows the 'true' surface density of quasars estimated in the various 
joint apparent magnitude and redshift bins using the maximum (a,b) and median (c,d) 
completeness estimators. The top panel for each estimated data set shows the predicted 
surface density using our best-fit LF model and the bottom panel that predicted by model 
B from Boyle et al. (1990), normalised appropriately. 
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Fig. 2. — This shows the quasar LF based on our best-fit model (solid line), extrapolated to 
z — 0, and the LF for Seyfert galaxies (Cheng et al. 1986) (*). The dotted line indicates 
the extrapolation of the quasar LF to magnitudes fainter than the limiting magnitude of 
the survey. The dashed line is model B from Boyle et al. (1990), extrapolated to z — 0, 
and the dotted-dashed line is the LF for soft X-ray selected AGN (Franceschini et al. 1994), 
extrapolated to z — 0. 



- 21 - 



Table 1. Data sets analysed and results. 





Data set 1 


Data set 2 


Data set 3 


Data set 1: subset 1 


Data set 1: subset 2 












> 0.7) 


(M B 


< -24.5) 


Detection technique 


UVX (0):163 


UVX (-0.3):135 


Colour:17 


UVX (0):125 


UVX (-0.35):98 


UVX (0):60 


UVX (-0.35):55 




Var.:145 


AQD:49 


X-ray: 19 


Var.:117 


Var.:117 


Var.:55 


Var.:55 


M[l] 


154 


92 


18 


122.5 


108.5 


57.5 


55 


M[2] 


16 


49 


1 


8.5 


18.5 


4.5 


5 


No. of quasars 


170 


141 


19 


131 


127 


62 


60 


Trial functions 
















F(v) = 
















K(l - av),v < a^ 1 


258 (66%) 


295 (48%) 


29 (66%) 


198 (63%) 


197 (50%) 


94 (64 %) 


91 (60%) 


0,v > a- 1 
















K(l - v) a 


190 (89%) 


302 (47%) 


21 (90%) 


141 (89%) 


154 (64%) 


67 (90%) 


66 (83%) 


K exp(— av 2 ) 


271 (63%) 


198 (71%) 


31 (61%) 


214 (58%) 


219 (45%) 


100 (60%) 


96 (57%) 


K8(v - a) 


172 (99%) 


197 (72%) 


19 (100%) 


132 (95%) 


131 (75%) 


62 (97%) 


61 (90%) 


K,v < a 


171 a (99%) 


262 (54%) 


18 a ( ■■■ ) 


176 (71%) 


175 (56%) 


60 a (100%) 


59 a (93%) 


0, v > a 
















Lower limit 


172 (99%) 


197 (72%) 


19 (100 %) 


132 (95%) 


131 (75%) 


62 (97%) 


61 (90%) 



a Solution has a > 1; this is the value for J F(y)dv (i.e. the upper limit is 1 and not a). 
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Table 2. 'Best-fit' parameters for evolutionary models derived using maximum and 

median completeness estimators. 



Luminosity Evolution a (3 M* kr> P(x 2 ) 



model 




Max 


Med 


Max 


Med 


Max 


Med 


Max 


Med 


Max 


Med 


Max 


Med 


2 power law 


(l + z) fci 


-3.50 


-3.51 


-1.24 


-1.33 


-23.38 


-23.60 


2.29 


2.20 






0.8732 


0.7087 


2 power law 


e k L r 


-3.56 


-3.56 


-1.26 


-1.35 


-22.26 


-22.51 


4.32 


4.12 






0.8572 


0.6721 


1 power law 


(1 + z) k L 


-1.83 


-1.83 






-22.71 


-22.99 


2.93 


2.73 






0.0048 


0.0033 


1 power law 


e k L T 


-1.85 


-1.84 






-16.62 


-21.39 


5.15 


4.72 






0.0078 


0.0040 


2 power law 


(1 + z) k t 


-3.35 


-3.15 


-1.00 


-0.86 


-22.51 


-21.97 


2.95 


3.36 


-1.06 


-1.79 


0.9385 


0.9064 


2 power law 


e k L T 


-3.42 


-3.33 


-1.07 


-1.00 


-21.23 


-20.75 


5.39 


5.92 


-1.08 


-1.72 


0.9392 


0.9391 


Boyle model B 


{l + z) k L 


-3.87 


-3.87 


-1.32 


-1.32 


-22.37 


-22.37 


3.20 


3.20 






0.0410 


0.0041 


Boyle model M 


{l + z) k L 


-3.90 


-3.90 


-1.34 


-1.34 


-22.59 


-22.59 


3.01 


3.01 


0.69 


0.69 


0.0162 


0.0008 


Boyle model N 


{l + z) k L 


-3.89 


-3.89 


-1.29 


-1.29 


-22.79 


-22.79 


2.78 


2.78 


1.70 


1.70 


0.0006 


4 x 10~ 6 



